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ABSTRACT 

We present a direct determination of the stellar metallicity in the close pair galaxy NGC 4038 (D 
= 20 Mpc) based on the quantitative analysis of moderate resolution KMOS/VLT spectra of three 
super star clusters (SSCs). The method adopted in our analysis has been developed and optimised 
to measure accurate metallicities from atomic lines in the J-band of single red supergiant (RSG) or 
RSG-dominated star clusters. Hence, our metallicity measurements are not affected by the biases and 
poorly understood systematics inherent to strong line H II methods which are routinely applied to 
massive data sets of galaxies. We find [Z]= +0.07 ± 0.03 and compare our measurements to H II 
strong line calibrations. Our abundances and literature data suggest the presence of a flat metallicity 
gradient, which can be explained as redistribution of metal-rich gas following the strong interaction. 
Subject headings: galaxies: individual: NGC 4038/39 — galaxies: abundances — galaxies: star clus¬ 
ters: general — methods: analytical — techniques: spectroscopic 


1. INTRODUCTION 

The metallicity of a galaxy is moderated by the cy¬ 
cling of chemically processed material by stars and any 
gas exchange between the galaxy and the environment. 
The central metallicity is correlated with a galaxy mass, a 
relation which holds information about galaxy formation 
and evolution ([Lequeux et al.|1979 Tremonti et al.|2004 


Maiolino et al. 12008 


2012). Additionally, 


Kudritzki et al. 
the variation ol the metallicity of the galaxy with the dis¬ 
tance from the centre keeps track of the complex dynam- 


processes affect metallicity gradients (e.g. 

Searle 19711 

Zaritsky et al.| 

19941 IGarnett et al.l|1997t (Frantzos & 

Boissier|20001 |C 

hiappini et al. 120011 |Fu et al. 120091 IPilk- 

ington et ai.|20H 

1 Mott et al.|20131 Kudritzki et al.|2015p. 

Obtaining reliat 

)le metallicities in galaxies is crucial to 


measure and interpret the behaviour of metallicity with 
radial distance and the mass-metallicity relation. Unfor¬ 
tunately, robust metallicity measurements in galaxies are 
notoriously problematic to obtain. Metallicity of star- 
burst or star-forming galaxies is routinely measured from 
H II region emission lines and two main methods are em¬ 
ployed: the Tg-based method and the strong line method. 
The Tg-based method uses the flux ratio of auroral to 
strong lines of the same ion to measure the electron tem- 

[T9M| 


perature of the gas dRubin et al. 
Stasihskal 120051 Andrews &: Martini 


Lee et al. 2004 

M| r 


However, 
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temperature sensitive lines are often too weak to be de¬ 
tected in faint distant galaxies and their measurement is 


ically in the metal-rich regime i 

Stasihska||2005 

Bresolin 

et al. 20051 

lErcolano et al. 12010 

|Zurita & Breso 

lin|2012[ 

Cazak et a 
not be dete 

. 12015 
rminec 

). When the e 

. one has to ref 

ectron temperature can- 
3ort to abundance indica- 


tors based on more readily observable lines. Such strong 
line methods are based on the ratio of the fluxes of the 


stronges t forbidden lines of typically Q and H (Pagel et| 
al.|1979 Skillman]1989 McGaugh|1994 ) which are more 
easily detected than the weak auroral lines across a wide 
range of metallicity. Unfortunately, well-known but un¬ 
explained systematic uncertainties, that can amount up 
to ~0.7 dex, plague the determination of the metallicity 
of extragalactic H II regions from nebular sp ectroscopy 
(Kewley Sz Ellison||2008 iBresolin et al. 2009). As a con- 
sequence, the chernical abundances derived from strong¬ 
line methods display large systematic differences when 
applied to the same observational data; and only rela¬ 
tive metallicity comparisons appear to be reliable if the 
same calibration is used. 

Very promising tracers of the present-day abundances 
in star-forming galaxies are evolved massive stars. In¬ 
deed, a number of nearby galaxies have been studied us¬ 
ing blue supergiants (BSGs) and results indicate excel¬ 
lent agreement with abundances obtained from the Tg- 


based method in H II regions (see Kudritzki et al. 2012 
20131120141 iHosek et al. 20141 and references therein). 


Tv growing body of evidence indicates that accurate 
metallicities over large distance scales as for BSGs can 
be derived also from red supergiant stars (RSGs). RSGs 
are young ( < 20 Myr) and extremely brig ht stars, i.e. 
10^-10® Lq (Humphreys & Davidson 1979). Their flux 
peaks at ~ l/i m, therefore they are among the most lu¬ 
minous objects in a galaxy in the near-IR and are ideal 
candidates for directly measuring extragalactic abun¬ 
dances. RSGs are also very cool stars, with tempera¬ 
tures ranging from 3000 to 5000 K, hence their spectr a 
show numerous absorption features (Allard et al.||200'0). 
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Previous techniques to measure metallicity from RSG 
spectra concentrated on the TJ-band, where high reso¬ 
lution (R[A/(5A]~ 20000) observations are needed to iso¬ 
late diagnostic atomic lines from molecular absorption 
(e.g. |Davies et al.||2009|). Even with the largest available 
telescopes, the need tor high resolution translates pro¬ 
hibitively large exposure times for individual objects at 
a distance >1 Mpc. In contrast, the required observing 
time can be significantly reduced if one focuses on a nar¬ 
row region in the J-band, where the dominant spectral 
features are isolated atomic lines of Fe and the a ele¬ 
ments (Ti, Si, and Mg). Indeed, in this spectral window 
accurate abundances can be measured even at moderate 
resolution (R[A/5A]~ 3000). 


The J-ban d method was initially introduced by Davies 
et al. (2010) for individual RSGs in the M ilky Way and 


has been extended and largely tested by 


( 2014a I in the association Perseus OB-1. 

2015p checked the validity of this method at lower metal- 


Gazak et al. 

Davies et al. 


licit y in the Magellanic G louds using the VLT/XShooter 


and Patrick et al. 
tion/analysis met 


(2015) accurately tailored the reduc- 
10 d tor KMOS observations of RSGs 


in NGC 6822. Gazak et al. (2015) obtained metallic- 


ities from RSGs across the disk ot NGC 300, a spiral 
galaxy beyond the Local Group, finding a striking agree¬ 
ment with the metallicities recovered from BSG stars and 
H Il-region auroral line measurements. 

Interestingly, this technique can be also applied to un¬ 
resolved star cluste rs rather than individual stars (Gazak 
et al. 2013 2014b). In merging and starburst galaxies 
intense star formation activity triggers the formation of 
super star clusters (SSCs), agglomerates of millions of 
young (< 50 Myr) stars. SSCs have masses > lO^M^r^ 
and are extremely compact (with radii < 5 pc; Portegies 
Zwart et al. 2010). Once a SSC reaches an age of ~ 7 
Myr, the most massive stars which have not yet exploded 
as supernovae will be in the RSG phase. For a cluster 
with an initial mass of 10® Mq, there may be more than a 
hundred RSGs present which dominate the cluster’s light 
output in the near-IR, contributing 90% - 95% of the of 
the near-IR flux. As their spectra are all very similar in 
the effective temperature (Tes) range around 4000 K, the 
combined spectrum can be analyse d in the same way as a 
single RSG spectrum, as shown by Gazak et al. (2014b). 
Therefore, for SSCs older than 7 Myr the J-band tech- 
nique can be used to measure metallicity at far greater 
distances than is possible for single super giants. 


Follow ing this line of investigation, Gazak et al. 


analysed two young SSCs in M 83 ( at 4.5 Mpc, 


.'him et a l. j 200 3 ) and NGC 6946 (at 5.9 Mpc, Karachent- 
sev et al. 2000), finding metallicities ~1.5-2.0 and ~U.5 
X solar, respectively. This paper further extends the 
observational baseline and presents quantitative J-band 
spectroscopy of three SSCs in the close pair galaxy 
NGC 4838. 

NGC 4038 is the main component of the Antennae 
system (NGC 4038/39), the closest (~20 Mpc)[^ and 
youngest example of an ongoing major merger, invmving 
two gas-rich disk galaxies that began to collide ~ 200- 


®|Riess et al.|||201l| estimated a distance modulus to the Anten¬ 
nae gaIaxTes''orrirM=31. 66 ± 0.08 (~ 22.3 Mpc) from optical and 
infrared observations of Cepheid variables with the the Wide Field 
Camera 3 (WFC3) on the Hubble Space Telescope (HST). 


400 Myr ago (I Barnes 1988 Mihos et al. 1993). Galaxy 
mergers, and their resulting starbursts, are one of the ba¬ 
sic b uilding blocks of struct ure formation in the universe 
(e.g.,|Baron & White||1987|) and represent an ideal labo- 
ratories tor close up studies of the physical processes that 
were important at the peak of cosm ic star formation (de 
Ravel et al. 2009 Bundy et al. 2009). As such, the Anten- 
nae has, over the years, been the favourite target for sev¬ 
eral multi wavelength studies of the effect of tidal interac¬ 
tion (e.g., Whitmore et al.|1999 Ehbbianoetalj20M Hi- 


bbard ^ al.||2005[ [Gilbert & Graham||2007[ [Brandi et al. 

20091 IWhitinore et al.|2010[ |Klaas et al.|2010| [Whitmore 


et al. [20T4) and numerous N-body and hydrodynamical 


simulations (e.g., 

Toomre & Toomrell 19721 

Barnes 

1988 

Mihos et al.||1993 

'I'eyssier et al.|2U10 Karl 

et al.|2010). 

'Lhe bodies ot 

formation (~20 IV 
an IR luminosity 

he galax 

I© yr-'; 

af log Lie 

ies are sites ot e: 

Zhang et al. 11200 
= 10.76, which 

ctensive star 

.) producing 
s is an order 


of magnitude lower than ultraluminous infrared galax¬ 
ies, but still a factor o f ^^5 higher than noni nteracting 
galaxy pairs (see, e.g., jKennicutt et al.||1987|. Most of 
the star formation in this colliding galaxy pair occurs in 
the form of SSCs dWhitmore & Schweizer||1995l iMirabell 


SSCs dWhitmore fc Schweizer||i995| iMirabe 
Whitmore°°er°ar 19991 [Wilson et al.||2000| 


et al. 1998 ^__ __^ 

with masses up to a few times lO^Mfr^ which are dis¬ 


tributed throughout the galaxy (Zhang et al. 2001). 

Observational studies demonstrate that m interacting 
galaxies, the strong metallicity gradient observed in iso¬ 
lated spirals, can be disrupted by gas flows of metal- 
poor gas from the outer regions t owards the centre of 
the galaxy (e.g. Kewley et al. 2010). 

Here we analyse KMOS/VLi' spectra of three SSCs to 
measure the central metallicity and the metallicity gradi¬ 
ent across the disk of NGC 4038. This paper is organised 
in the following way. In Section!^ we summarise the ob¬ 
servations and data reduction. We outline our analysis 
procedure in Section]^ Our results are presented in Sec¬ 
tion U] We summarise our main conclusions in Section [sl 

2. TARGET SELECTION, OBSERVATIONS AND DATA 
REDUCTION 

Our sample consists of three RSG dominated SSCs, 
whose coordinates, luminosities, colours, ages are listed 
in Table[^together with other useful information. Target 
SSCs with the appropriate luminosity and colour s were 
selected from opt ical and near-IR photometry by Whit-| 


more et al. (2010) from the Advanced Camera for Surveys 
(AGS) ana the Near Infrared Camera and Multi-Object 
Spectrometer (NICMOS) mounted on the HST. Their 
spatial location across the galaxy is shown in Figure 
The observations we re carried out with KMOS/VLT 
(Sharpies et al.||2013) in April 2014 in visitor mode (PI: 
Kudritzki: 093.B-0023), with a total exposure time of 
6000s split in 20 subsequent exposures. 

KMOS is a spectrograph equipped with 24 deployable 
integral-field units (IFUs) that can be allocated within 
a 7.2' diameter field-of-view (FoV). Each IFU covers a 
projected area on the sky of about 2.8"x 2.8", which 
is sampled by an array of 14x14 spatial pixels (spaxels) 
each with an angular size of 0.2". The 24 IFUs are man¬ 
aged by three identical spectrographs, each one handling 
eight IFUs. The observations were performed in nod-to- 
sky mode with the YJ grating, covering the 1.00-1.35 fim 
spectral range with a nominal resolution of R~3600 at 
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Fig. 1.— The targeted SSCs are overlayed over a HST-ACS F814W image of NGC 4038 retrieved from the Hubble Legacy Archive 
(http://hla.stscl.edu/hlaview.html). 


TABLE 1 _ 

SSCs TARGETED. DATA ARE FROM [WhITMORE ET AL.|([2010[| 


ID 

R.A. 

Decl. 

Mv 

mjA) 

E{B - V) 

Logr/yr 

Mass (Mq) 

35897 

12:01:50.4453 

-18:52:14.223 

-11.6 

15.8 

0.00 

7.6 

4.5 xlO® 

36731 

12:01:55.9896 

-18:52:12.985 

-12.5 

17.3 

0.04 

7.6 

1.1 xlO® 

50776 

12:01:51.3963 

-18:51:47.562 

-14.4 

16.4 

0.04 

6.8 

1.1 xlO® 


J magnitudes are from aperture photometry on archival Wide Field Camera 3 images. The 
typical uncertainty on the J magnitude is ~ 0.1 mags. 


the band centre. Observations were carried out using the 
standard AB AB-like object-sky sequence (i.e. one sky 
frame for each object frame) in which we offset by 5" to 
sky, and each observation was dithered by up to 0.2". 
During the observations, the average J-band seeing was 
approximately 1.0". In addition to science observations, 
a standard set of calibration frames were obtained. Since 
we require high precision absorption line spectroscopy, 
we observed a telluric standard with the arms in the sci¬ 
ence configuration (i.e. using the observational template 
that allows users to observe a standard star in each IFU 
allocated to a science target). During the observations 
7 SSCs were observed, b ut only three of tho se have suf¬ 
ficient SNR (> 100, see|Gazak et al.||2014a| to be used 
in our analysis. To reduce the data, we used the stan¬ 
dard recipes provided by the by the Software Pac kage for 
Astronomic al Reduction with KMOS (SPARK; 


et al. 20131. KMOS IFU data cubes were flat fielded. 


Davies 


wavelength calibrated, and telluric corrected using the 


standard KMOS/esorex routines. 

As can be seen in Figure there are significant varia¬ 
tions in both spectral resolution and wavelength calibra¬ 
tion across the FoV of each IFU as measured from sky 
emission lines (see below for more details). Left uncor¬ 
rected, these can introduce sky and telluric cancellation 
errors into the final spectrum, which can be the source of 
substantial problems for precision absorption line spec¬ 
troscopy such as that presented in this paper. 

We correct for th ese effects with a pr ocess we call, kmo- 
genization (see also Gazak et al. ]20l^ . We first take the 
rectified science and sky images prior to sky subtraction. 
At each spaxel in the IFU we fit Gaussian profiles to the 
sky lines and we use them to determine both the spec¬ 
tral resolution and the higher order wavelength calibra¬ 
tion as a function of spatial position. We then smooth 
the spectra at each position on the IFU down to a lower 
resolution, set to be R = 320C[3 We then extract the 

® Though there are regions on the IFU where the spectral reso- 
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Fig. 2. — Variation of spectral resolution and wavelength cali¬ 
bration across the spatial pixels of a typical KMOS IFU measured 
from night sky lines. The top panel presents the reconstructed im¬ 
age obtained in a 300 sec exposure. The middle and bottom panels 
show the variation of spectral resolution and the wavelength shifts 
(in pixel) across the spaxels, respectively. 


lution is below R = 3200, the pixels within the extraction aperture 
are always above this value. 


spectrum of the science target from an aperture with 
radius 1.5 spaxels around the flux peak. This narrow 
aperture minimises errors due to spatial non-uniformity, 
at the slight expense of discarding flux from the wings of 
the point-spread function. When the source is extracted, 
the wavelength axis is updated to include the higher or¬ 
der correction determined from the sky lines without in¬ 
terpolating the spectrum onto a new wavelength axis, as 
this can introduce numerical noise. The associated sky 
spectrum is extracted in the same way, and is subtracted 
from the science spectrum. 

Since the wavelength calibration can change slightly 
throughout the night due to the rotation of the instru¬ 
ment, each sky-subtracted science spectrum can have a 
different wavelength axis. To account for this, the nor¬ 
malised spectra are co-added onto a master wavelength 
axis whereby the flux at wavelength in the master 
spectrum is a weighted mean of the fluxes in the indi¬ 
vidual spectra which have wavelengths between Ai_i and 
Ai+i. Rejection of outliers is also performed at this stage, 
to eliminate cosmic rays and bad pixels. 

Finally, we experimented with altering the output 
spectral resolution to test the robustness of our method. 
We found that the best-fitting parameters (see next sec¬ 
tion) were stable to well within the fitting errors as long 
as the output resolution was greater than R = 3100. 


3. ANALYSIS 


Metallicity [Z] (normalised to Solar values, [Z] = log 
(Z/ZqI^ effective temperature (Teff), surface gravity 
(log g) and microturbulence (^) have been derived as ex¬ 


(e-S- 

Davies et 

al. 20101 

Gazak et al. 

2014b|a 

Davies 

et 

al. 

21)15 Patrick et al. 

20151 Gazak et al. UM 

tlTTR 

particular, the studies ot 

Davies et al. (I2U1UD, 

Gazak et 

al. 

(|2014a 

1, and 

Gazak et al. (2U15p demonstra) 

te the ap- 


chemical enrichment. 

Atmospheric parameters and metallicity were derived 
by comparing the observed spectra with a grid of single¬ 
star synthetic spectra degraded to the same spectral res¬ 
olution as those observed. Mod el atmospheres were cal - 
culated with the MARCS code ( Gustafson et al.||2008[) 


where the range of parameters are dehned in Table [2 
The synthetic spectra were computed usi ng the update 
vers ion o f the SIU code, as described in [Bergemann et| 
SI (2012[). Departures from local thermodynamic equi¬ 


librium (LTE) for Fe I, Mg I, Si I, and Ti I lines were 


also included (Bergemann et al. 2012[ 2013 2015). All 


other lines including the weak molecular contributions 
are calculated in LTE. The best ht model has been de¬ 
rived through a -minimisation between the observed 
spectrum and a template spectrum at each point in the 
model grid, taking into account possible shifts and vari¬ 
ations in instrumental spectral resolution between the 
observed spectra and the models. The methodology for 
finding the best-fitting model is described in detail by 


We measure Fe, Mg, Si, and Ti abundances from individual 
lines and assume this is representative of the metallicity Z. While 
the assumption of a solar-scaled composition appears reasonable, 
we note that the errors deriving from the assumption of a solar- 
scaled composition rather than an [o/FeJ-enhanced one has a little 
impact on the metallicity measurements, i.e. well within the quoted 
errors on metallicity. 
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TABLE 2 TABLE 3 

Model grid Spectral Fits 


Parameter 

Units 

Grid Min 

Grid Max 

Grid Step 

SSG 

Teff(K) 

log g (dex) 

? (km/s) 

Z (dex) 

Tefl 

K 

3400 

4400 

100 

35897 

3890 ± 50 

0.4 ± 0.2 

1.7 ± 0.2 

0.01 ± 0.07 

log g 

dex 

-1.0 

± 1.0 

0.5 

36731 

3750 ± 50 

0.6 ± 0.2 

2.0 ± 0.2 

0.09 ± 0.08 

[Z] 

dex 

-1.00 

±1.00 

0.25 

50776 

3770 ± 50 

0.5 ± 0.2 

2.0 ± 0.2 

0.11 ± 0.07 


km 

1.0 

5.0 

1.0 







Gazak et al. (2014a) and Gazak et al. (12015 ) which anal- 
yse 11 RSGs in Perseus UB-1 ([ZJ = -0.04 ± 0.08) and 
27 RGGs in NGC 300 ([Z] = -0.03 ± 0.05), respectively. 
Briefly, we measured best fit parameters by isolating 2D 
planes in each parameter pair combination which contain 
the minimum. As a result, two parameters are locked to 
the best fit values in each of the six slices, while allowing 
the other two parameters to vary. We thus constructed a 
2D plane of minimum values as a function of the two 
fixed parameters and draw contours of equal y^. This 
procedure is repeated for all the planes, and the mea¬ 
surements of each parameter are finally averaged to de¬ 
rive best fit parameters. An example of this is shown 
in Figure]^ Figure]^ shows the degeneracy between [Z] 
and log g, while demonstrating ho w soundly con s trained 
are Tpff and £. We ref er again to |Gazak et al.| (|2014a|) 
and |Gazak et al.| (|2015|) for a thorough discussiori of the 
sensitivity of diagnostic lines to the free parameters (e.g. 
Tgff, log g, ■Cj ^-nd [Z]) and the complete error analysis. 
The best-fitting parameters are listed in Table along 
with their associated uncertainties. 

4. RESULTS 
4.1. Metallicity 

Figure [^displays fits of the metal lines with the best 
fitting model spectra in our diagnostic spectral window. 
The match to the stellar lines used to derive metallicity 
(labelled in the same Figure) is very good, and also the 
fit to the unresolved features of the pseudo-continuum is 
remarkable. The spectral feature at 1.205/im is due to a 
poorly removed telluric line. 

Table summarises the derived stellar parameters and 
metallicities. The average metallicity for our sample of 
3 SSGs is [Z] = -1-0.07 ± 0.03 (cr=0.05 dex), in agree¬ 
ment with previous measurements based on diffe rent di¬ 
agnos t ics. S olar metallicity has been inferred by |MengeT| 
et al. ( 2002 ) from the analysis of the metallicity sensi¬ 
tive Mg I l ine at 8806.8A in a handful of SSGs. 


tian et al. 


Bas- 


measure the strong emission lines in 
H II region j^^]~to derive metallicities of 16 SSGs, find¬ 
ing a slightly super-solar value (see Figure E- They 
use the equivalent width ratios of the collisional excited 
[O II]A3727 and [O III]A4959, 5007 lines relative to the 
H/? Balmer series recombination lines to estimate the gas- 
phase oxygen abundance, i.e. the R23 lin e ratio. More¬ 
over, they ad opt the analysis method from [Kobulnicky 
Kewley|2004| (hereafter KK04) and a solar abund ances of 
12 - 1 - log(U/H )0 = 8.89 ^Edmunds fc Pagd 1984), which 
is 0.2 higher than ou r adopted solar metall ici ty, i.e. 12 - 1 - 
log(O/H )0 = 8.69 ( Asplund et al.[ 2009). Bastian et 


As the lifetimes of RSGs are <50 Myr, their metallicities are 

ex£ected_t^j3e_^e2resentative_^^^e_^j£udB_^|om_wMchMhe^_are 

tormed. 


iVlOS/VLT spectra. However, they do not have available 
the [O II]A3727 emission line str ength to measure the 


0/H ratio. Therefore, they use the Vacca & Gonti (19921 
calibration that relies on the R3 ratio. Since this cali- 
bration does not include the contribution from another 
ionised speci es of oxygen ( i e. _ [O II]), the abundances 


presented by Bastian et al. (2006) should be taken with 
caution, as stated by the authors themselves. In the fol¬ 
lowing, we do not discuss those abundances furtheip^ 

4.2. Metallicity gradient 

Any sizeable galaxy-galaxy interaction affects the 
metallicity distribution of the galaxies involved. Obser¬ 
vational studies of interacting or close pair galaxies have 
shown that these undergo nuclear metal dilution owing 


gas-phase metallicity gradients dLee et al.|2004 

Tranche 

et al. 

120071 IChien et al.|20071 

Kewley et al.||2UUbt|Kupke 

et al. 

l20081|Ellison et al.|20081 

Michel-Dansac et al. 

2008 

Kewlev et al. 12010 Rich et al. 

2012 Sanchez et al. 

2014 

Rosa et al.||2014|). 


analysing the influence of different levels of interactions 
in the metallicity distribution of galaxies. [Kewley et al.| 
( 2010 ) present the the first systematic analysis of metal- 
hcity gradients in close pairs. They obtain spectra of star 
forming regions in eight galaxies in close pair systems and 
find that the metallicity gradients are signific antly shal¬ 
low er tha n the gradients in isolated spirals. [Krabbe ^ 
^ (|2008|) study the kinematics and physical properties 
of tne minor merger AM2306-721. They report a clear 
metallicity gradient across the disc of the main galaxy, 
while the secondary, less massive companion shows a rel¬ 
atively homogeneous oxygen abundance. A nearly flat 
radial gradient was measu red for both compon e nts of the 


system A M 2322-821 by (Krabbe et al. 2011). Werk et 
al. ( 2011 ) find that the interacting galaxies m their spec¬ 


troscopic sample have flat oxygen abundance gradients 
out to large projected radii. A flat metallicity gradient 
has also been reported for NGG 1512, a barred spiral 
in a clo se interaction with a blue compact dwarf com¬ 
panion ( jBresolin et al. 2012 ) and NGC 92, the larges t 
galaxy in Robert 's quartet (Torres-Flores et al. 2014). 
Rosa et al. (2014) measure oxygen abundance gradients 
from H II regions located in 11 galaxies in eight systems 
of close pairs from Gemini/GMOS spectra, finding metal- 


Note that two of the SSGs reanaivsed in|Bastian et al.|||2006 


- - - _ , flip-- -, 

were also observed in (Bastian e t al.l (|2009||. and one o n our targei 
(36731) is in common with the |Bastian et ah] (|2006|l sa mple; i.e. 
their compl ex 3. For the SSC in common with our sample, (Bastian] 
|et al.|||2006|l estimate a metallicity of Z=0. 45 Zr?,, with an estimated 
intrinsic uncertainty of ±0.2 in log(0/H) lEdmunds &: Page lll984(l 
and a solar oxygen abundance of 12± log (U/tlj0= S.yi jMejierl 
|1985[(, which is compatible with our measurement. 
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Fig. 3.— Spectra of the targeted SSCs in the J-band window, along with the best-fit spectra (red lines). 



UK) 




Fig. 4.— Planes through the model grid showing the degeneracy 
between the parameters for the analysed SSCs . The black contours 
show contours of equal and are drawn at levels of Xmin~^ 

3, 5, 10). The red contour is drawn at =Xmin~^^ which indicates 
our la uncertainty. The x-axis limits on the right-hand plot are 
the minimum and maximum possible values of log g allowed by the 
object temperature and luminosity. 


licity gradients significantly flatter than those observed 
in typ ical isolated spiral galaxies. Finally, [Sanchez et al.| 
(2014) derive radial gradients of the oxygen abundance 


in ionised gas in 306 nearby galaxies observed by the 
CALIFA survey. Using a large and homogeneous sam¬ 
ple of more than 40 mergers/interacting systems, they 
find statistical evidence of a flattening in the abundance 
gradients in interacting systems at any interaction stage. 

The absence of an abundance gradient in interacting 
galaxies is explained by invoking efficient mixing of low- 
metallicity gas from the outer parts with the metal-rich 
gas of the centre of the galaxy. S moothed particle hy¬ 
drody nami c merger simulations of [Barnes &: Hernquist] 
(1996) and Mihos & Hernquist ( 1996| ) predict that in in¬ 
teracting galaxies the strong tidal interaction during en¬ 
counters can aid disk to develop bars. In such barred 
disks, low-metallicity gas from the outskirts can effi¬ 
ciently flow towards the central regions at higher metal- 
licity funnelled by bar inst abilities, flattening the initial 
radial metallicity gradient (Hibbard fc: van Gorkom|1996[ 


Georgakakis et al. |2000[ Kampazzo et al. 


ipke et 


.||2(I05l llono et al. 1120051 [Emonts et al.||20()6| |Martin 

i06| ICullen et al.||2007^ 


al 

200 ^^ 
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Fig. 5.— Direct metallicities from the SSC analysis (large 
sq uares) are shown togeth er with the H II region metal abundances 
by [Bastian et al.| ( |2009| l as a function of galactocentric radius. 
Small symbols are tne a bundances obtained from d ifferent strong¬ 
line calibratio ns to the [Bastian et al.|[2009| data: |Kobulnick^& 
Kewl e vl2004l (KK04: red emp ty circles, adopted bylBa stian et al. 
2UU!j rTRewley fc Dopita|2002|(KD02: blue filled circles I, jlVic Oiaugn 
lyyi (jviyi; teal empty triangles), and|Pettini & Pagel|zUU4|(FFU4; 
green filled triangles). 


We plot our direct metallicity measurements against 
the radial distance from the centre in Figure Galac¬ 
tocentric radii are computed using the direct distance 
to the galaxy centre, where the galaxy centre is defined 
using Two Micron All Sky Survey (2MASS) images. Dis¬ 
tances take inclination into account (f=56.56° ) using the 
optical diameters from Lauberts & Valentijn| (|1989 ). A 
weighted linear regression to our SSG data yields a cen¬ 
tral metallicity of +0.17 ± 0.23 dex and a fiat gradient 
(-0.02 ± 0.03 dex kpc“^), consistent with observations 
and N-body simulations of interacting systems. 


4.3. Strong-line Abundances 

Even though the original observations do not include 
all the different metallicity indicators, we can still com- 
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pare our direct SSC r netallicities to strong-line calibra¬ 
tions, as applied to the|Bastian et al. (2009) datrp^ SSCs 
rnetallicities refer to the combined abundances~Mg, Si, 
Ti, and Fe (see Section where the individual abun¬ 
dances are scaled relative to the solar abundance pat¬ 
tern). In comparing the rnetallicities from the SSCs with 
the oxygen abundances from the H II regions, we assume 
that the oxygen abundances of the SSCs scale with metal¬ 
licity, and that the solar metallicity value corr esponds to 
12 -I- logfO/H )0 = 8.69 (jAsplund et al.||2009|). 

Figure [^illustrates how tour dittererir stroiigd ine cali¬ 
brations, as applied to the Bastian et al. (2009) sample, 
compare to our metallicitiesl We consider the R23 = 
([O II]A3727 -I- [O III]AA495 9,5007)/H/3 indic ator, with 
the theoretical predictions by Mc Caugh|1991 (M91) an d 
KK04, i.e. the calibration used by Bastian et al. (2009), 
the calibrati on for the [N II]A6583/|U il|A3727 diagnostic 
presented in Kewley & Dopita 2002 (KD02), and 03N2 
= log ([O IIIJ A5007/H/t) / ([IN IIJ A6584/Hq;) , empirically 
calibrated on H II regi ons with Tg-based rnetallicities by 
Pettini & Pagel||2004| (PP04). The estimated accuracy 
of the all these calibrations is ~0.10 - 0.15 dex. We use 
the coeffic i ents o f the polynomial from Table 3 in |Kewley| 
'& Ellison (2008) to convert rnetallicities from the KK04 
calibration into a metallicity that is consistent with other 
calibrations using a third-order polynomial. 

The abundances that result from the application of 
these different indicators and calibrations, together with 
the direct SSC abundances derived in Section [^ are 
shown as a function of the distance from the centre of 
NCC 4038 in Figure A linear least-square fit to the 
data is also shown for each method. While all the consid¬ 
ered methods imply a flat metallicity gradient, both the 
KK04 and KD02 calibrations provide rnetallicities that 
are larger than those from SSCs, while the PP04 cali¬ 
bration give rnetallicities that are lower than ours. Fig¬ 
ure [5] also shows that the M91-based rnetallicities better 
agree to the direct metallicity determination. However, 
we stress that only the KK04 and M91 calibration s rely 
on the set of em ission lines actually observed in |Bas-| 
tian et al.| (|2009|) data, i.e. [O II1A3727, [O III1A4959, 
[U 1IIJA5UU7, H/t; while both the PP04 and KD02 cali¬ 
brations include diagnosti c lines (i.e. [N II]A6 583, Ha) 
not present in the original Bastian et al. (2009|) spectra. 
Moreover, we note that a!ll the abundaiices from Fig- 
urej^agree within ±2cr. As a result, we caution that the 
comparison presented in Figure is purely qualitative, 
because we are not comparing actual measurements but 
using instead empirically derived metallicity conversions 
with their associated uncertainties. 

The data suggest the presence of a flat abundance gra¬ 
dient, however they sample only a small range in dis¬ 
tances from the galaxy centre. Indeed, the number of 
SSCs observed should be increased considerably, espe¬ 
cially in the inner regions, in order to test for radial vari¬ 
ations and confirm whether the metallicity gradient is as 


shallow as found by Bastian et al. (2009) and suggested 
by our data. 

Knowledge of the chemical composition of galaxies is 
fundamental to trace back the history of cosmic chemi¬ 
cal enrichment and understand the processes at work in 
galaxy formation and evolution. Chemical abundances 
in extragalactic environments are commonly based on 
H II region optical emission-line ratios. Numerous rela¬ 
tions have been proposed to convert diagnost ic emission 
line ratios in to metallicity (see for a review Kewley & 
Dopita 2002). Nonetheless, comparisons between abun¬ 
dances obtained from different calibrations show a sys¬ 
tematic offset in metallicity estimates, that can amoun t 
up to ~0.7 dex ( Bresolin||2008 Kewley fc Ellison||2008| ). 

We have developed a iiew method to avoid these cali- 
bration issues by performing quantitative spectroscopy of 
RSC stars or RSC-dominated SSCs in external galaxies. 
In this paper we use this J-band technique on KMOS 
spectra of three SSCs in NGC 4038 to measure precise 
rnetal licities (see Daviesetalj201^ Gazak et al. 2014b|a 


2015| Davies et al. 2015” Ratrick et al. 2015() We hire 
an average metallicity of [ZJ = -|-0.07=[: 0.03 ((t= 0.05 
dex). Given the uncertainties, this is in g ood agreement 
with t he results f rom H II regions data by [Bastian et al. 


(2009), when the McCaugh (1991) calibration is used to 
determine abundances, furthermore, we find no evidence 
for a metallicity gradient. However, a larger systematic 
study of RSGs is needed to assess the presence (or not) 
of a metallicity gradient among the young population 
within this merger galaxy pair. 

With the multi-object near IR spectrographs such as 
KMOS/VLT and MOSFIRE/Keck we can now investi- 


from individual RSC stars (i.e. 

Davies et al. 112010 

2015 

Patrick et al. 

2015 

Gazak et aF 

2015). Using the 

same 


as 0.10 dex out to the astonishing distance of ~20Mpc 
in less than one night of observations. This opens new 
windows for extragalactic spectroscopy. Indeed, the J- 
band method will allow us to quantitatively study the 
chemical evolution of galaxies -up to the Co ma cluster - 
in a w ay similar to current Galactic studies (Evans et al. 
20111 , when the next generation of extremely large tele- 


scopes, equipped with adaptive optics supported near IR 
multi-object spectrographs, will be available to the com¬ 
munity. 

Based on observations made with ESO Telescopes at 
the La Silla Paranal Observatory under programme ID 
093.B-0023. We are very grateful to the anonymous ref¬ 
eree who helped us improve our manuscript. RPK and 
JZG gratefully acknowledge support by the National Sci¬ 
ence foundation (NSF) under grant AST-1108906. BP 
thanks the CNRS Programme National de Physique Stel- 
laire for financial support. 

Facilities: VLT(KMOS). 
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